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ABSTRACT 

We show that the observed mass-to-light ratio of galaxy clusters increases with 
cluster temperature as expected from cosmological simulations. Contrary to previous 
observational suggestions, we find a mild but robust increase of M/L from poor (T 
~ 1 — 2 keV) to rich (T ~ 12 keV) clusters; over this range, the mean M/Ly increases 
by a factor of ~ 2. The best fit relation satisfies M/L v = (170 ± 30)T^ 0A h at z = 0, 
with a large scatter. This trend confirms predictions from cosmological simulations 
which show that the richest clusters are antibiased, with a higher ratio of mass per unit 
light than average. The antibias increases with cluster temperature. The effect is caused 
by the relatively older age of the high-density clusters, where light has declined more 
significantly than average since their earlier formation time. Combining the current 
observations with simulations, we find a global value of M/Ly = 240 ± 50/i, and a 
corresponding mass density of the universe of O m = 0.17 ± 0.05. 

Subject headings: galaxies: clusters of - cosmology: theory - cosmology: observations 
- dark matter - large-scale structure of the universe 



Introduction 



The mass-to-light ratio of clusters of galaxies has been used for decades in the classical M/L 
method of estimating the mass density of the universe. In this method, the average ratio of the 
observed mass to light of clusters of galaxies, which represent systems of relatively large scales (~ 1 
Mpc), is assumed to be a fair representation of the mean M/L of the universe. This ratio is then 
multiplied by the observed luminosity density of the universe to yield the universal mass density. 
When applied to rich clusters, with average M/Lb — 300/i, the total mass density of the universe 
amounts to £l m ~ 0.2 (where £l m is the mass density in units of the critical density; Zwicky 1957; 
Abell 1965; Ostriker, Peebles, & Yahil 1974; Bahcall 1977; Faber & Gallagher 1979; Trimble 1987; 
Peebles 1993; Bahcall, Lubin, & Dorman 1995; Carlberg et al. 1996; Carlberg, Yee, & Ellingson 
1997, and references therein). 
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A fundamental assumption in this determination, however, is that the mass-to-light ratio of 
clusters is a fair representation of the universal value. If M/L of clusters is larger or smaller than 
the universal mean, the resulting Q m will be an over- or under- estimate, respectively. In general, 
if mass follows light on large scales (M/L ~ constant), the galaxy distribution is considered to 
be unbiased with respect to mass. If mass is distributed more diffusely than light, as is generally 
believed, then the galaxy distribution is biased (i.e., more clustered than mass), and the above O m 
determination would be an underestimate. 

Observations of galaxies, groups, and clusters of galaxies show that the mass-to-light function 
increases with scale up to a few hundred Kpc (Rubin & Ford 1970; Roberts & Rots 1973; Ostriker, 
Peebles, & Yahil 1974; Einasto, Kassik, & Saar 1974; Davis et al. 1980; Zaritsky et al. 1993; Fischer 
et al. 2000), but then flattens on larger scales (Bahcall, Lubin, & Dorman 1995; Fischer et al. 
2000). A comparison with high resolution cosmological simulations (Bahcall et al. 2000) reveals 
an excellent agreement with the observed M/L function: the simulated M/L function increases 
with scale on small scales (< 0.5 Mpc) and flattens to a mean constant value on large scales, as 
suggested by the data. 

However, the simulations show that while M/L flattens on average on large scales, high over- 
density regions, such as rich clusters of galaxies, exhibit higher M/L values than average, while 
low density regions have lower M/L ratios (at a given scale); high density regions are antibiased in 
M/L, with mass more strongly concentrated than light (a higher M/L ratio) than average. Bahcall 
et al. (2000) show that the M/L antibias exists in simulated clusters of galaxies: higher density 
clusters (richer, higher mass and temperature clusters) should exhibit larger M/L ratios (especially 
in the blue and visual luminosity bands) than poorer, lower mass and temperature clusters. The 
expected effect is approximately a factor of two increase from poor to rich clusters. It is caused, 
as seen in the hydrodynamic simulations, by the age of the systems: high density clusters form at 
earlier times than low density clusters; their luminosities have thus declined more significantly by 
the present time (if no significant recent star formation occurred), yielding a larger M/L ratio, on 
average, for the richer clusters. 

Observationally, however, cluster M/L ratios have been suggested to show no dependence on 
cluster overdensity (as represented by cluster richness, mass, or temperature) (e.g., Carlberg, Yee, & 
Ellingson 1997; Hradecky et al. 2000). In this paper, we investigate the observed M/L dependence 
on cluster overdensity (temperature) by using a well-studied sample of clusters that range from 
poor to rich systems. We find that M/L increases, on average, with cluster temperature, in good 
agreement with the prediction of cosmological simulations. 

2. Dependence on Cluster Temperature 

We investigate the dependence of M/L on cluster overdensity by using a reliable sample of 
clusters that ranges from poor to rich systems. The cluster overdensity is represented by the 
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observed temperature of the cluster, T; the overdensity within a given radius is proportional to T: 
(Ap/p) c ; oc (M/R 3 ) c i oc T. As revealed by cosmological simulations (Bahcall et al. 2000), the M/L 
ratio of clusters is expected to increase with T by a factor of ~ 2 (in M/Lb and M/Ly) from poor 
(T ~ 1 - 2 keV) to rich (T ~ 10 keV) clusters. 

Current observations of galaxy clusters generally claim no observed dependence of M/L on 
cluster temperature or other cluster parameters (e.g., Carlberg, Yee, & Ellingson 1997; Hradecky 
et al. 2000). However, the expected effect is not large, and the observational uncertainties in M/L 
are significant. It is important therefore to use a reliable and consistent sample of clusters that 
spans a wide range of cluster temperatures. 

We select clusters with measured temperatures in the range of ~ 1 to 12 keV. All clusters 
have measured M/L ratios observed within radii of typically 0.5 to 1.5 h^ 1 Mpc. The cluster M/L 
measurements are based on a few self-consistent subsamples and methods. We analyze different 
subsamples separately in order to test the reliability of the results. The sample includes clusters 
with mass determinations from gravitational lensing observations (Fischer Sz Tyson 1997; Hoekstra, 
Franx, & Kuijken 2000; Squires et al. 1996a, 1996b, 1997; Tyson & Fischer 1995) and X-ray 
observations (from a single, carefully executed method; Hradecky et al. 2000). A sample of velocity 
dispersion virial mass determinations from a single method (Carlberg, Yee, & Ellingson 1997) is 
also shown for comparison, although this method differs somewhat in that it determines cluster 
virial masses, within the larger (and less certain) virial radii of ~ 1 — 1.5 h^ 1 Mpc. We repeat our 
analysis both with and without this virial subsample and find consistent results. We also repeat 
the analysis for individual cluster subsamples, e.g., lensing versus X-ray cluster subsamples and 
high versus low redshift cluster subsamples; all samples yield consistent results. 

In addition to the cluster sample, we also include the recently observed ensemble average M/L 
ratio of 50 groups of galaxies determined from weak gravitational lensing by Hoekstra et al. (2001). 
We convert the mean velocity dispersion for this group sample, 274+sg km/s, to a mean temperature 
using the observed velocity-temperature relation for clusters and groups: V r = (332±52)T^ e y km/s 
(Lubin & Bahcall 1993). This yields a mean group temperature of T = 0.73 keV, consistent with 
observations of typical groups (~ 1 keV). The total sample we use, 21 systems, is listed in Table 
1. Also listed are the redshifts, the observed temperatures and M/L ratios, the luminosity band 
of the observations (mostly V band), the relevant radius R for the observed parameters, and the 
relevant references. A Hubble constant of Hq = lOO/i km/s/Mpc is used throughout. 

Two corrections are applied to the sample to bring all measurements to the same consistent 
system. First, the sample contains clusters at redshifts from z ~ 0.02 to 0.8. Since cluster lumi- 
nosities evolve with redshift, we correct all luminosities to z = 0. Given the observed evolution 
in Ly (Kelson et al. 1997), Ly ~ iq(0-3z±w%) ^ a rj M/Ly values were corrected upward by a fac- 
tor of lo(°- 3z±10% ) to yield M/L v (z = 0). Similarly, M/L B values were corrected by a factor of 
10 (o.4z±io%) ( van Dokkum et al. 1998) and M/L r values were corrected by a factor of lo(°- 15z±10% ) 
(Carlberg et al. 1996). 
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Second, four of the 21 clusters are observed in the blue band luminosity, Lb, and six arc 
observed in the r band luminosity, L r (all others are observed in Ly). We convert Lb to Ly 
at z = using Ly ~ 1.15 L#, based on the typical cluster color of B — V ~ 0.8 (and solar 
(f? — V)q = 0.65). Similarly, we convert L r to Ly using Ly ~ 0.935 L r , following (5 — r) c ; ~ 0.915 
and {B — r) & = 0.692 (Jorgensen, Franx, & Kjaergaard 1995). These corrections have no significant 
impact on the final results. In addition, all cluster values (and the evolution correction) are corrected 
to the same cosmology of qo = 0.15 (following van Dokkum et al. 1998). The resulting M/Ly values 
at z = are listed in Table 1. The error bars represent the convolved error bars of the observations 
and the evolutionary correction. 

The sample includes cluster M/L ratios measured within radii of ~ 0.5 to 1.5 hr x Mpc. The 
results are not significantly affected by the different radii; we test this by comparing a subsample 
of clusters measured at a fixed radius (0.5 h~ l Mpc) and find consistent results. The cluster M/L 
ratio is nearly constant over this range of radii within a given cluster. 

The observed M/Ly(z = 0) ratios are plotted as a function of cluster temperature for all 
clusters in Figure 1. While the scatter is large, a clear correlation of increasing M/Ly with T is 
apparent. The best fit linear and power-law relations are: 



These best fit relations are shown in Figure 1. Also shown in the figure is the relation predicted 
by the cosmological simulations (Bahcall et al. 2000). The agreement between data and simulations 
is excellent; both show an increase in M/Ly by a factor of ~ 2 on average, from ~ 200/i to 450/i, 
as T increases from ~ 1 to 12 keV. While the simulations M / L(T) trend is insensitive to Q m , 
the normalization, or absolute value of M/L at a given T, is proportional to VL m . The simulation 
results plotted in Figure 1, which best fit the data, correspond to Cl m = 0.17 (using the observed 
(M/Ly) critical — 1400/i as the simulation normalization; see Bahcall et al. 2000). 

Figure 2 represents the same data as Figure 1, but with the data binned in several temperature 
bins (3 keV each). The average M/Ly for each bin is shown. The best fit relations, also plotted, 
are consistent with those obtained for the full sample (Figure 1, eq. [1], eq. [2]). The best fit results 
and their x 2 are summarized in Table 2. 

To test the sensitivity of the results to the cluster selection method and to the radius used, we 
repeat the analysis for two subsamples of the full sample: 1) we omit the clusters with dynamical 
determination of virial masses (to R v ~ 1 — 1.5 h~ l Mpc) based on velocity dispersion data 
(Carlberg, Yee, & Ellingson 1997); and 2 ) we use only clusters measured within a fixed radius of 
0.5 hT l Mpc. The best fit results obtained for each of these subsamples are consistent with the 




z=0 



= 142 ± 32 + (23 ± 5) T keV h = 142 ± 32 [1 + (0.16 ± 0.05) T keV ] h 



(1) 




(2) 
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results of the full sample; they are summarized in Table 2. The results for the R = 0.5 h Mpc 
subsample are shown in Figure 3. 

We further test the results by analyzing other subsamples of clusters: clusters with gravitational 
lensing masses versus clusters with X-ray mass determinations, as well as clusters at low redshift 
(z < 0.1) versus clusters at intermediate redshift (z ~ 0.15 to 0.33). All subsamples yield consistent 
results (within la). These tests reinforce the robustness of the results. 

The above results are in excellent agreement with the prediction of cosmological simulations, 
showing an increase in cluster M/Ly with cluster temperature (or overdensity) . 

3. Summary and Conclusions 

We show that the mass-to-light ratio of clusters of galaxies increases, on average, with cluster 
temperature. Contrary to previous observational suggestions of a constant M/L ratio for clusters - 
independent of temperature, velocity dispersion, or other cluster parameters - we show a consistent 
trend of increasing M/Ly with temperature from groups and poor clusters (at T ~ 1 keV) to the 
richest clusters (T ~ 12 keV); over this range, the mean M/Ly ratio increases by a factor of ~ 2. 
The best fit relations are given by equation (1) and equation (2) (and Table 2, Figures 1-3). 

This trend is in excellent agreement with the prediction of cosmological simulations (Bahcall 
et al. 2000) (see Figures 1 - 3 for comparisons between data and simulations). While the scatter is 
large, the underlying trend is clear. 

The observed correlation is important for several reasons. 

1) It confirms the predictions originally made from cosmological simulations (Bahcall et al. 
2000) of increasing cluster M/L ratio with cluster overdensity (temperature). As seen in the 
simulations, this trend is caused by the "age effect": higher density systems formed earlier and 
their luminosities, especially in the blue and visual bands, have decreased more significantly by the 
present time relative to the younger, lower density systems. This results in larger M/L ratios for 
the older, higher temperature clusters. 

2) The simulated results, confirmed by the current observations, show that rich clusters are 
antibiased in their M/L ratios - that is, rich clusters have higher M/L ratios than average. This 
implies that mass is more concentrated than light, on average, in rich clusters (Bahcall et al. 2000; 
see also Jing, Mo, Sz Borner 1998). This antibias increases with cluster richness (temperature). 

3) Antibias in the cluster M/L ratio directly affects the measurement of the mass density 
of the universe, £l m . Classically, £l m determinations use the observed M/L ratio of the richest 
clusters (which are easiest to observe) and assume it to be representative of the global value. 
The observed antibias, however, shows that this method overestimates Q m , since rich clusters 
have larger M/L ratios than average. Comparing the observed dependence of M/L on T with 
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results from the cosmological simulations (Bahcall et al. 2000), shown in Figures 1-3 [using 
{M / Ly) critical = 1400/i ± 20% for the simulation normalization (see Bahcall et al. 2000)], we find 
a global mass density parameter of £l m = 0.17 ± 0.05. The mean representative mass-to- light ratio 
for the universe is M/Ly = 240 ± 50ft, comparable to that exhibited by groups and poor clusters. 

REFERENCES 

Abell, G. O. 1965, ARA&A, 3, 1 
Bahcall, N. A. 1977, ARA&A, 15, 505 

Bahcall, N. A., Lubin, L. M., & Dorman, V. 1995, ApJ, 447, L81 

Bahcall, N. A., Cen, R., Dave, R., Ostriker, J. P., & Yu, Q. 2000, ApJ, 541, 1 

Carlberg, R. G., Yee, H. K., Ellingson, E., Abraham, R., Gravel, P., Morris, S., & Pritchet, C. J. 
1996, ApJ, 462, 32 

Carlberg, R. G., Yee, H. K., & Ellingson, E. 1997, ApJ, 478, 462 

Davis, M., Tonry, J., Huchra, J., & Latham, D. 1980, ApJ, 238, L113 

Donahue, M., Voit, G. M., Gioia, I., Luppino, G., Hughes, J. P., & Stocke, J. T. 1998, ApJ, 502, 
550 

Einasto, J., Kaasik, A., & Saar, E. 1974, Nature, 250, 309 
Faber, S. M., & Gallagher, J. S. 1979, ARA&A, 17, 135 
Fischer, P., & Tyson, J. A. 1997, AJ, 114, 14 
Fischer, P. et al. 2000, AJ, 120, 1198 

Hjorth, J., Oukbir, J., & van Kampen, E. 1998, MNRAS, 298, LI 
Hoekstra, H., Franx, M., & Kuijken, K. 2000, ApJ, 532, 88 
Hoekstra, H. et al. 2001, ApJ, 548, L5 

Hradecky, V., Jones, C., Donnelly, R. H., Djorgovski, S. G., Gal, R. R., & Odewahn, S. C. 2000, 
ApJ, 543, 521 

Irwin, J. A., & Bregman, J. N. 2000, ApJ, 538, 543 

Jing, Y. P., Mo, H. J., & Borner, G. 1998, ApJ, 494, 1 

Jorgensen, L, Franx, M., & Kjaergaard, P. 1995, MNRAS, 273, 1097 



- 7- 



Kelson, D. D., van Dokkum, P. G., Franx, M. , Illingworth, G. D., & Fabricant, D. 1997, ApJ, 478, 
L13 

Lubin, L. M., & Bahcall, N. A. 1993, ApJ, 415, L17 
Mushotzky, R. F., & Scharf, C. A. 1997, ApJ, 482, L13 
Ostriker, J. P., Peebles, P. J. E., & Yahil, A. 1974, ApJ, 193, LI 

Peebles, P. J. E. 1993, Principles of Physical Cosmology (Princeton: Princeton Univ. Press) 
Roberts, M. S., & Rots, A. H. 1973, ApJ, 186, L95 
Rubin, V. C., & Ford, W. K. 1970, ApJ, 159, 379 

Squires, G., Kaiser, N., Babul, A., Fahlman, G., Woods, D., Neumann, D. M., k, Bohringer, H. 
1996a, ApJ, 461, 572 

Squires, G., Kaiser, N., Fahlman, G., Babul, A., k Woods, D. 1996b, ApJ, 469, 73 

Squires, G., Neumann, D. M., Kaiser, N., Arnaud, M., Babul, A., Bohringer, H., Fahlman, G., k 
Woods, D. 1997, ApJ, 482, 648 

Trimble, V. 1987, ARA&A, 25, 425 

Tyson, J. A., & Fischer, P. 1995, ApJ, 446, L55 

van Dokkum, P., Franx, M, Kelson, D., k Illingworth, G. 1998, ApJ, 504, L17 
Zaritsky, D., Smith, R., Frenk, C., k White, S. D. M. 1993, ApJ, 405, 464 
Zwicky, F. 1957, Morphological Astronomy (Berlin: Springer) 



This preprint was prepared with the AAS IATjrjX macros v5.0. 



-8- 



Table 1. Cluster parameters. 



Object 


z 


T 


M/L 


Band 


R 


M/L v 


Ref. 






(kcV) 


(h Mq/Lq) 




(h- 1 Mpc) 


(z = 0) 




A262 


0.016 


2.3 ±0.2 


180 ± 64 


V 


0.5 


180 ± 64 


7 


A426 


0.018 


6.2 ±0.4 


154 ± 68 


V 


0.5 


154 ± 68 


7 


A478 


0.09 


8.4t°-j 


276l 2 ?| 


V 


0.5 


288l^ 7 


7 


A1795 


0.062 


7.8 ± 1.0 


468lJ« 


V 


0.5 


479li|| 


7 


A2052 


0.035 


2.8 ±0.2 


198 ± 54 


V 


0.5 


201 ± 55 


7 


A2063 


0.035 


2.3 ±0.2 


190 ± 46 


V 


0.5 


193 ± 47 


7 


A2199 


0.030 


4.8 ±0.1 


294 ± 110 


V 


0.5 


297 ± 111 


7 


MKW4s 


0.028 


1.8 ±0.3 


202 ± 86 


V 


0.5 


204 ± 87 


7 


A2163 


0.20 


11.0 ±0.6 


300 ± 100 


V 


0.5 


331 ± 110 


8, 12 


MS 1054 


0.83 


12.3+22 


248 ± 34 


B 


0.5 


464 ± 74 


2, 5 


A2218 


0.18 


7.48toi? 


440 ± 80 


B 


0.4 


478 ± 87 


4, 10 


A1689 


0.18 


9-0218:18 


400 ± 60 


V 


1 


435 ± 65 


4, 13 


A2390 


0.23 


8.9018:?^ 


320 ± 90 


V 


1 


360 ± 101 


4, 11 


RXJ1347 


0.45 


1, O-+1.10 

-0.92 


200 ± 50 


B 


1 


264 ± 67 


4, 3 


MS0016 


0.55 


8.0 ± 1.0 


213 ± 74 


r 


1.32 


275 ± 96 


4, 1 


MS1358 


0.33 


6 50+ ' 68 
D - OU -0.64 


188 ± 40 


r 


1.15 


225 ± 48 


4, 1 


MS0440 


0.20 


c 0+1.3 
o - o -0.9 


314 ± 105 


r 


0.87 


360 ± 120 


9, 1 


MS0451 


0.54 


io.i7ii : y 


314 ± 105 


r 


1.45 


405 ± 136 


9, 1 


MS0839 


0.19 


4 in+0.36 
^' ±a -0.33 


317 ± 132 


r 


1.12 


362 ± 151 


9, 1 


MS1008 


0.31 


7.29+»-« 


247 ± 79 


r 


1.36 


294 ± 94 


9, 1 


50 groups 


0.33 


n 7 o+0.07 


191 ±83 


B 


0.5 to 1 


225 ± 98 


6 



References. — (1) Carlberg, Yee, & Ellingson 1997; (2) Donahue et al. 1998; (3) Fischer & 
Tyson 1997; (4) Hjorth, Oukbir, & van Kampcn 1998; (5) Hoekstra, Franx, & Kuijkcn 2000; (6) 
Hoekstra ct al. 2001; (7) Hradecky et al. 2000; (8) Irwin & Brcgman 2000; (9) Mushotzky & 
Scharf 1997; (10) Squires ct al. 1996a; (11) Squires ct al. 1996b; (12) Squires ct al. 1997; (13) 
Tyson & Fischer 1995. 
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Table 2. Best-fit M/L v - T relations. 



Sample 


M/L v = A T keV + B 


X 2 /dof 


M/L V = CT« v 


X 2 /^/ 


Full sample (21) 


(23 ± 5)T + (142 ± 32) 


0.8 


(173 ± 29)T°- 30±0 - OS 


1.0 


Subsample 1*(15) 


(21 ± 6)T + (153 ± 39) 


1.2 


(171±33)T - 30±01 ° 


1.3 


Subsamplc 2*(11) 


(26 ± 7)T + (129 ± 40) 


1.0 


(133 ± 35)T°- 44±0 - 14 


1.1 



* Subsample 1: excludes clusters with dynamical determination of virial masses 
*Subsamplc 2: clusters observed at a fixed radius (R = 0.5 h^ 1 Mpc) 
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Fig. 1. — Observed M/Ly (corrected to z = 0) vs. T for clusters of galaxies. The best fit linear and 
power law relations (eq. [1] and eq. [2]) are presented by the solid and dotted lines, respectively. 
The prediction from simulations is shown by the dashed line; its normalization corresponds to 

n m = o.i7. 
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Fig. 2.- 



Same as Figure 1 but for the data binned in 3 keV temperature bins. 



